THE MUTUAL INCLUSION IN A NONLOCAL 
COMPETITIVE LOTKA VOLTERRA SYSTEM 
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Abstract. We investigate the traveling front solutions of a nonlocal 
r—{ , Lotka Volterra system to illustrate the outcome of the competition 

between two species. The existence of the front solution is obtained 
through a new monotone iteration scheme, the uniqueness of the front 
solution corresponding to each propagation speed is proved by sliding 
domain method adapted to nonlocal systems, and the asymptotic decay 
rate of the fronts with critical and noncritical wave speeds is derived 
by a new method, which is different from the single equation case. The 
results demonstrate that in the long run, two weakly competing species 
can co-exist. 



1. INTRODUCTION 

In this paper we study the properties of traveling front solutions of the 
following nonlocal Lotka- Volterra competition system: 

(i.i) ( 5 = f*?-5+% 1 r fi ^° 1 ^ (M)elxl + 

v ' I vt = J * v — v + rv{\ — (12U — v), 

where '*' denotes convolution: J*w = J^J(x—y)w(t,y)dy. The nonnegative 
functions u(x,t) and v(x,t) are population densities of the two competing 
species. The constant r > is the relative growth rate of species v, and a±, 
a>2 > are interactive constants. The integration kernel J is a probability 
density function satisfying JL J = 1, being radial and nonnegative with finite 
support of nonzero measure. Depending on the conditions of a\ and a 2 , sys- 
tem (jl.ip may have four constant equilibrium states O(0, 0), A{1, 0), B(0, 1) 
and ; 1— a°a )• ^ e can then classify their stability as follows: 

(i) If < ai,a 2 < 1, CEO has four equilibria with di~°* , ) being 
asymptotically stable, and (0,0), (1,0), (0,1) being unstable; 

(ii) If < ai < 1 < ci2, (jl.ip has three equilibria: (1,0) is asymptotically 
stable, while (0,1) and (0,0) are unstable; 

(iii) If ax, a 2 > 1, then (1,0) and (0, 1) are both stable and ( ^""^ , x j_~" 2 2 ), 
(0, 0) are unstable; 
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(iv) If < 02 < 1 < a\, then (0,1) is asymptotically stable, while (0,0) 
and (1,0) are unstable. 

We are interested in the case (i) so the condition < 01,03 < 1 is as- 
sumed throughout the paper. The cases (ii) and (iv) are currently being 
investigated in [IB] . 

System (jl.ip describes the long range competition between two species, so 
a natural question is who will be the winner in the long run. We will address 
this question, and others such as the uniqueness, the monotonicity as well 
as the asymptotics by studying the traveling front solutions of (jl.ip . The 
traveling front solution of (jl.ip is a pair of nonnegative smooth functions 
(u(x,t),v(x,t)) = (u(x + ct),v(x + ct)) := (-u(£), u(£)), which satisfies the 
following boundary value problem 



where £ = x + ct, x € M, t £ M + is the moving coordinate and c > is the 
speed of the front. 

The nonlocal reaction diffusion equations and systems can be found in 
many applications such as cell biology, phase transition, ecology, and the 
neurons and neuronal network [H [2l El [HI [2Ql H31 [251 (ZD El L231 LT] - One 
of the important features of those models is the appearance of traveling 
front solutions, which is frequently used to study the transition between two 
equilibrium states, see ([H]). In recent years, there have been fruitful results 
concerning traveling front solutions of the nonlocal equations and systems. 
For example in [11] the authors systematically investigated the traveling 
front solutions of integro-differential equation 

(1.3) u t (x,t) = J*u(x,t) -u(x,t) + f(u(x,t)), (x,t)eRxR + 

and established the existence and uniqueness of the traveling front solution 
as well as the asymptotic behaviors for the ignition type, the bistable type 
and the KPP (Kolmogorov-Petrovsky-Piskounov) type nonlinearities. And 
in [9], the nonlocal diffusion-reaction equations of the form 



where Q is a bounded interval, was considered. In particular, the maximum 
principle as well as the sliding domain method were derived, and successfully 
applied to obtain the uniqueness and monotone behavior of positive solution 
of (jl.4p . Recently, Pan and Li et al [22] established the existence of the trav- 
eling wave fronts in nonlocal delayed reaction-diffusion systems by using the 
upper-lower solution method and Schauder's fixed point theorem. G. Lv [19] 
dealt with the asymptotic behavior of the traveling fronts and entire solu- 
tions for a nonlocal monostable equation. By using lower and upper travel- 
ing wave solutions, the authors in [25^ proved the existence of traveling wave 




(1.2) 




(1.4) 
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solutions for the two-species integro-differential Lotka-Volterra competition 
model in the form similar to For the recent work concerning the non- 

local reaction-diffusion equations and systems, refer to [25], ?, Ell [22l [24] 
and the references therein. 

We note that other methods, such as the homotopy continuation method 
[U [2], and the contraction mapping method [20], could also be used to 
establish the existence of traveling front solutions for nonlocal equations 
and systems. Here we would like to propose a new and easy monotone 
iteration scheme to derive the existence of traveling wave solutions of 
To this end, we first change (jl.2p into a local monotone type by setting 
u = 1 — u, v = v, so this leads to 

{J * u — u — cv! + {u — l){u — a\v) = 0, 
J * V — V — CV' + rv(l — (12 + CL2U — v) = 0, 
(u,v)(-oo) = (0,0), (u,v)(+oo) = (aik 2 ,k 2 ). 

Observing that if u is replaced by a±v in the second equation of (|1.5p . we 
recover the nonlocal KPP equation. And the existence, uniqueness as well as 
the asymptotic behaviors of the traveling front solutions are known recently 
([HI [HI H]). The construction of the upper and lower solutions for (jl.5p is 
based on this simple observation. 

The asymptotic behavior closely ties to the other properties of the front 
solutions such as their uniqueness and the stability. The asymptotics of the 
traveling wave solutions is investigated by several steps of integral inequal- 
ities. Our method differs from that in [U 119] . where Ikehara's Tauberian 
Theorem can be applied directly. In the current system the method fails due 
to the nontrivial coupling of the two components. We solve this problem 
by first studying the canonical form of the linearized system, then switching 
back to the original system to have the desired results. 

For the proof of the uniqueness of traveling front solutions of (ll.5p . we 
adapt the generalized sliding domain method |16j . Combining the sliding 
domain method and the nonlocal comparison principle we are able to es- 
tablish the uniqueness (up to a shift of the origin) of the traveling wave 
solutions for each front speed. 

The outline of the paper is as follows. In Section [2] we set up the upper 
and lower solutions and then apply monotone iteration scheme to derive the 
existence of the front solutions. In Section [3] we prove the strict monotonic- 
ity of the traveling front solutions by comparison principle and obtain the 
asymptotic decay rates, and derive the uniqueness of the front solutions by 
a generalized sliding domain method. 

2. EXISTENCE OF FRONT SOLUTIONS 

In this section, we first introduce some lemmas to aid set up the upper 
and lower solution pairs for (11.5p . The existence of traveling front solutions 
is then established by monotone iteration ([22]). For the rest of the paper 
the inequality between two vectors is component- wise. 
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Definition 1. A smooth function , £ G R is an upper solutions 

of (II. 5p if it satisfies 



(2.1) 



J * u — u — cu' + (u — l)(u — a\v) < 0, 
J * v — i? — cv' + rv(l — a 2 + — v ) <■ 



and the boundary conditions 

(2.2) (u,u)(-oo) > (0,0), (u,v)(+oo) > (a 1 k 2 ,k 2 ). 

A lower solution of (jl.5p is defined in a similar way by reversing the inequal- 
ities in (HJ]) and (j2~2"1) . 

The construction of upper and lower solutions pairs is based on the solu- 
tions of the nonlocal KPP equation of the form: 



(2.3) 



J * u — u — cu' + f(u) = 0, 
u(— oo) = 0, ii(+oo) = b > 0, 



where the smooth function / satisfies /(0) = f(b) = 0, /(s) > for s € 
(0, 6), /'(0) > and / is non-increasing near 6, the integral kernel J satisfies 
the conditions specified in section [TJ 
Let the constant c* be defined as 

c* =: min {a > : \ \ [ J&e^dx + /'(0) - ll j, 

then c* > and is finite ([23J) . We recall the following result: 

Lemma 2. Let f be defined as above, then for any c > c* system (|2.3p 
has a unique (up to a shift of origin) smooth and monotonically increasing 
solution w{£), £ € R. 

Proof. See [23], Theorem 2.1]. 

□ 

To construct the upper solution for the system (jl.5p . we begin with the 
following version of nonlocal KPP equation: 



(2.4) 



J * v — v — cv' + r(l — a2)v{l — jj-) = 0, 
v(— oo) = 0, u(+oo) = &2 



where corresponding to (|2.3p /(S) = rv{l — a2+a\a2V — v) >0for ; uG (0, /c2), 
recalling that k 2 = x ^~ a ^ ■ It is easy to see that /(0) = /(fc 2 ) = 0, /'(0) = 
r(l — a 2 ) > and f'(k 2 ) = — r(l — a 2 ) < 0. According to Lemma [21 for each 
fixed c > c* with 

(2.5) c* =: min f A > : - [ / J(£)e A? cfe + r(l - a 2 ) - ll ) > 0. 

A L7r -I J 

system (|2.4p has a unique (up to a translation of the origin) traveling front 
solution v satisfying the given boundary conditions. 
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Define 

then we have the following result, 

Lemma 3. Let c* > be defined as in (|2.5|) then for each fixed c > c* 
is a smooth upper solution for system (II. 5h . 



Proof. On the boundary one has 

u(— oo) \ _ / \ / u(+oo) \ _ f CL\kl 



v( — oo) y V y ' \ U(+oo) J V 
As for the it component, we have 

J *u — u — cv! + (u — 1) (u — ai u) 
= ai( J * tJ — u — cv') + (aiW — l)(ai^ — a\v) 
= —a\rv(l — a2 + aia2V — v) 
= — a\r{X — a2)v{k2 — v)<0 

We then verify the second component in (jl.5j) . 

J * v — v — cv ' + rv(l — a2 + a\a2V — v) 
= —rv(l — a2 + aia2V — v)+ rv(l — 02 + a\a2V — v) 
= 

Thus (u,v) forms a smooth upper-solution for (|1.5[) . 

□ 

We next construct the lower solution pair for system (jl.5p . We choose a 
number Z according to one of the following conditions: 



(2.7) 



(2 



' 1). < I < max{^, 1+r[1 _ a2 l +k2[ai _ r) }i{0 < r < a x . 

2) .0<K TT ^ifa l <r<^- 

3) .0 < I < max{ , r ~ ai i\ . -rrrr T l if r > i 

' Lai(ra 2 -1)' l+r(l— a 2 ) J 02' 

and work with another nonlocal KPP equation: 

f J * « - u - cv' + r(l - 02)v(l - 1 7! 1 a a 2 2 ^ ) = 0, 
1 ^(- oo )= ^(+ oo ) = T i ^ 2 7<^. 2 

Corresponding to the notions in Lemma [21 

f(v) = r(l-a2)v(l- 1 —^v)>0 
V 1 — a 2 / 

fOT « e (0, j^j), and /(0) = /(t^) = 0, /'(0) = r(l - a 2 ) > and 
f'i lS^a 1 ) = ~ r (l ~~ °2) < 0. Noting we have the same c* as in the previous 
Lemma. 



THE COMPETITION INCLUSION 



6 



For each fixed c > c*, define 
where t>(£) is solution of (j2.8j) . 



Lemma 4. For eac/i /jxed c > c* , (j2.9j) is a smooth lower solution of system 

dm 



Proof. On the boundary one has 

TZ(— oo) 



and 



due to the fact 



u(+oo) 
u(+oo) 



v —oo 



1 — 0,2 



l—a±a2l I <- 
1-02 J — 



1— aia2l 



Furthermore, 



1 — Q2 
1— a\d2l 



< k 2 . 



J * v — v — cv_ + rv(l — a>2 + — 1>) 

/ 1 — a\02l \ > 
— r(l — a2j!M 1 2) + vyyi — ci2 + a\a2ly_ — vj 



1 — a 2 



0, 



Next, we check lower solution with respect to u, the computation as follows: 

J * u — u — cu_ + (u — 1) (u — a\v) 
= a\l(J * v — v — a/) + (a±lv — l)(a\ly_ — a\v) 



a\l 



r(l — e^JV I 1 : v 



+ aiv(a\lv — !)(/ — 1) 



1 — a 2 

= aiv{-rl{\ -a 2 ) -1 + 1 + [-rl(aia 2 l - 1) + aiZ(7 - l)]v} 
> 



by the condition (|2.7p . Hence, the conclusion of the lemma follows. □ 



We next show that the upper and lower solutions as obtained in Lemmas 
[3] and [5] are ordered. This can be achieved by shifting the upper solution 
to the left along the axis, but first a comparison of the asymptotic rates of 
the upper and lower solutions at infinities ia needed. Consider the following 
characteristic equations at ±co: 

M_(A) = f J(£)e Af d£ - cA + /'(0) - 1, at - oo, 

and 
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M + (A) = f J(£)e A€ d£ - cA + /' (b) - 1 at + oo. 

It is easy to verify thatM_(A) and M + (A) ([23] and [1]) have the following 
properties: 

(a) The characteristic equation M_ (A) for nonlocal reaction-diffusion equa- 
tions (12. 3|) (at — oo) has two positive roots Ai(c) < A2(c) for c > c* and one 
double root A* (c) > for c = c* and no real root for < c < c* ; 

(b) The characteristic equation M + (A) (at +oo) has one negative root 
As(c) and one positive root X^c) for any c > 0; 

The next lemma specifies the asymptotic behaviors of the front solutions 
of the KPP equations ([19] and [4]): 

Lemma 5. The traveling front solution w(£) as in Lemma^has the follow- 
ing asymptotic behaviors: 

For the critical front with speed c = c* , 

w (£) = b w £e x *^ + o(£e A *^), £ ^ -oo 

w (£) = b- d w e X3 ^ + o(e A3(c * )€ ), £ +oo. 
and for the noncritical front with c > c* , 

w {£) = a w e Xl ^ + o(e Al ^), £ -> -oo; 

w(£) = b- c w e X3{c) t + o(e A3(c) «), £ -> +oo; 
where a w , c w , d w are positive constants and b w is negative. 

By shifting the upper solution far enough to the left, the upper- and 
lower-solutions as derived in Lemma [3] and Lemma [Hare ordered. 

Lemma 6. Let c > c* be fixed and (u, v) T , (u,v) T be the upper and lower 
solutions defined in (|2.6|) and (|2.9p . then there exists a number r > such 
that 

( J)(e + r)>(|)(0 /or£ G M. 

Proof. We only show the lemma for c > c* since the one for c = c* is similar. 
By Lemma 

(2.10) (l)(0=( ai f i 1 ) eXl ' + O{eXl " ) 
and 

(2.11) ( I ) (0 = ( ) e A ^ + o(e A ^) 
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as £ — >• — oo; and 



(2.12) | ; «)= X - XX + 



and 

(2.13) [Z)i& = {7 Jr^j-^* j«~+*~> 

as £ — >• +oo, where Ai > is the smaller positive root of the characteristic 
equation 

A/_(A)= / J(Oe A? de-cA + r(l-a 2 )-l, 

JR 

and A3 < is the negative root of the characteristic equation 

M + (A)= / J(£)e A? d£ - cA - r(l - a 2 ) - 1, 
and Ai, A2, B±, B2 are positive constants. 

Since (|2.4|) is translation invariant, I> T (f;) =: v(£ + r) is also a solution of 
()2.4p for any f £l. It then follows that (v? ,v T ) T (£) is an upper solution for 
system (|1.5|) . For the asymptotic behavior of (u T ,v T ) T '(£) at —00, we can 
simply replace {a\A\, Ai) T by (aiAi, Ai) T e XlT in (|2.10p . Now we choose 
r > large enough such that 



n 



Ai J \ B x 

Then there exists a sufficiently large N\ > such that 

On the other hand, the boundary conditions of upper and lower solutions 
at +00 also imply that there exists a number N2 > such that 

(2.15) (rUMiD^ 61 ™ 

We next show that the inequalities (|2.14p and (|2.15p also hold on the 
interval [— iVijiV^]. There are two possible cases to deal with: 
Case 1. If we already have 

(2-16) (%%)>( 7$) on[- Nl ,N 2 ]. 



we then let r = t and have the conclusion. 

Case 2. There exists a point £0 £ ( — Ni, N2) such that 

tf'Ko) A „ / «(&) 



(2 ' 17> I J £ V ,(«o) 

with strict inequality holding for at least one of the two components. 
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In this case, we can use the sliding domain method [16]. We first shift 
(u T ,v T ) T (£) to the left by increasing r until we can find a, t\ > t > such 
that Cu Tl (b>I> Tl (£)) T > (u(OiH(0) T on the interval [-Ni, N 2 - (n - r)]. 
we then shift (u 71 ^), v T1 (£)) T back to the right by increasing n to some 
t 2 > t such that one of the branches of the upper solution tangents to its 
counterpart of the lower solution at some point £2 in the interval (— N\ + 
i~2, N 2 — (ti — r)). On the endpoints of the interval (— N\ +r 2 , ./V2 — (ri — r)), 
we still have (^(£),^ T2 (£)) T > Cu(£))2^(£)) T - In summary, we now have 
^ 2 (6) = u(6) and u T2 (£) > u(0,u r2 (0 > v(0 for £ G (-iVi + r 2 ,iV2 - 

(n-r)). 

Let := {if 2 ,v T2 ) T (£) - {u,v) T (t) and F = F 2 ) T = ((u-l)(u- 

a\v), rv(l — a 2 + 02^ — V )) T - F° r £ € (— iVj. + T2, ./V2 — (ri — r)) we have 

(2.18) 

J*^! -101 - Ctl^ + + Cl^l,«)wi + + (2W 2 )W 2 < 0, 

J * w 2 - w 2 - cw' 2 + ^(m+Cs^i,^! + ^(u, v + Uw 2 )w 2 < 

for some Q G [0, 1], i = 1,2,3,4. Since the above system is monotone we 
further have 

dFi _ 
(2.19) J * w\ — w\ — cw\ + — — (u + (iWi, v)w\ < 0. 

du 

That £2 is a global minimum point for u>i(£) implies w[(^ 2 ) = 0- Hence at 

( = ^2 we have 

(2.20) 

dFi 

J * hji(£ 2 ) - 'Wife) - cw[{i 2 ) + -^-{u + CiWi,v)w 1 (C2) = J * ^1(6) > 0. 

This contradiction shows that such £2 does not exist, and therefore we can 
further decrease t 2 to r. This shows that £0 does not exist either, we there- 
fore have 

SUMS!) 

□ 

We still use (U, v) T to denote the shifted upper solution as given in Lemma 
El With the ordered upper and lower solutions in Lemma O and the mono- 
tone iteration scheme |22j . we have the following result on existence of the 
traveling wave solution of (|1.5p . 



Theorem 7. Assume that < a\,a 2 < 1 and (12. Jjj , for each c > c* 
with c* being defined in (12.5p . system \T75\ exists a traveling wave solution 
(u(x,t),v(x,t)) = (u(x + ct),v(x + ct)) connecting (0,0) and {a\k 2 ,k 2 ). 

Proof. Applying the monotone iteration method |22] of the ordered upper 
and lower solution pairs obtained in Lemma [3] and Lemma [H we then con- 
clude the Theorem. □ 
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Remark 8. As pointed out in the introduction, there are several methods 
to show the existence of the traveling wave solutions, and the above pro- 
posed construction of the upper and lower solution is one of the easiest; and 
furthermore, the method can also provide some insight into the asymptotic 
behavior of the traveling front solutions, as it will be shown in the next 
section. 



3. ASYMPTOTIC BEHAVIOR AND UNIQUENESS OF THE FRONT SOLUTIONS 

The monotonicity of the front solutions is a direct consequence of the 
nonlocal comparison principle and the monotone iteration. 

Theorem 9. The traveling front solution of (|1.5j) derived in Theorem^ is 
strictly monotonically increasing from — oo to +oo. 

Proof. By the monotone iteration process [22], the traveling wave solution 
U(£) = (u(^),v(^)) T is increasing for £ € R, then its derivative satisfies 
(9i(0,92(0) T = U>(0>0, and 

(3.1) J * gi - g x - cg[ + (u, v)g\ + (u, v)g 2 = 0, 



dF 2 dF 2 

(3.2) J * 92 - 92 - cg' 2 + -g^i u ^ v )9i + ~q^^ v )d2 = 0, 

and 

(3.3) G?i(£),<? 2 (£)) T > 0, ( 5l , 52 ) T (±oo) = 0. 

The local monotone structure of (13.11) and (13.21) as well as the the maximum 



principle ([10]) implies that (<?i, g 2 ) T ^£) > for £ £ R. This concludes the 
strict monotonicity of the traveling wave solutions. □ 

We next derive the asymptotic behaviors of traveling front solution at 
±oo. Such result for the single nonlocal equation with KPP or bistable 
nonlinearty is already known ([H HH]), but for the system with nontrivial 
coupling such as (jl.ip it is new. 

Since the asymptotic behavior at — oo can be obtained by a straightfor- 
ward comparison, we only concentrate on deriving the asymptotic behavior 
at +oo. 

To get the asymptotic decay rate of the traveling wave solutions at £ — > 
+oo, we set c > c* and use the vector form to derive the results. Firstly, 
changing the system (11.5P by the transformations: 



u = a\k 2 — n, v = k 2 — v 
and reversing the sign of £, we will have the following system 

-J * u + u — cu' + (—ki — u){a\v — u) = 0, 
-J * v + v — cv' + r(k 2 — v)(—a 2 u + v) = 



(3.4) 
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and yet the boundary conditions are still 

(3.5) (l)(-co)-(S). (»)<+=.)-(■£ 

We only need to study the asymptotic behavior of the system (|3.4p and 
()3,5p at £ — > — oo. Once this is done, on switching back to the original 
variables we have the asymptotic behavior at +00. 

Lemma 10. Let c> c* be fixed then there exists a number 7 > such that 
the solutions of (|3.4p and (|3.5p satisfy 



u (0 \ - n,.--i\ 



, < O(e^) as £ — > -00. 
v{C) J 

Remark 11. This Lemma says that the traveling front solution of system 
(jl.5p approach the equilibrium (a\k2, ^2) at least exponentially fast for £ — > 
+00. 

Proof. We divide the proof into three steps: 
Step 1. is integrable for £ close to —00. 

For £ negatively large the system (I3.4p satisfies 



c[%) = -J*( ») + (») + { (~h - u)(a lV - u) 



(3.6) 

v' J u ^yvJ'\vJ' \ r(k2 — v)(—a2U + v) 

_j t[ u\(u\{ h-e - ai k 1 -e\fu 



v J V v J V ~ a 2 r k2 — e ^2 — e J \ v 
The above inequality is true because the limit of the Jacobian of the vector 
((—fei — u)(ai7; — u),r(k2 — v)(—a2U + v)) as £ — > —00 satisfies component- 
wisely the inequality 

k\ —a\k\ \ / k\ — t —a\k\ — e 
-a 2 r/c 2 r/c 2 / ~ \ _a 2 r ^2 — e rk 2 - e 
for a sufficiently small positive number e and sufficiently large £ < 0. 

Integrating the above inequality (I3.6P from y to £ (close to —00) we have 

$)-($))*jf(-"(S) + (5»* 



+ ( kl 7 e - a l kl ~ e ) f( u ~ ){s)ds 

x -a 2 rk 2 - e rk 2 - e J \ v 1 



or equivalently 
(3.7) 



ki — e 


—a\k\ — 


a 2 rk 2 - e 


rk2 — e 


m \ 


( 5(y) ^ 


m ) 





i? ( I ) w* 



£ «=(( ~S )-( Z! ))+£(•>*( 5 - S ))* 
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Since u,v are both bounded, the left hand side of the above inequality is 
bounded as long as the last integral in the right hand side is finite. 

By Fubini's Theorem and Lebesgue's Dominated Convergence Theorem 

%,(!)-(!))«„ = -jf^ w (( = ) (.,-( = ) (.-„),«. 

zJ(z) / ( ~, J (s — tz)dtdzds 



zj(z) f (~) (a-tz)dtdz 



as y — > — oo, therefore 



is bounded. 

This shows that 



J*\ I - ~ MM* 



k\ — e —aiki - e \ /" / n 
-a 2 rA:2 - e rk 2 - e J / V v 



(s)ds 



is bounded, which implies that 



(s)ds 



is also bounded. This is due to the fact that the matrix 

k\ — e —a\k\ - e \ 1 _ 1 / rk 2 — e a 2 rk 2 + e 
-a 2 rk 2 — e rk 2 - e J D \ a±ki + e fcj — e 

has all positive entries and D = {ki — e){rk 2 — e) — [— a\k\ — e] [— a 2 rk 2 — e] > 
for sufficiently small < e < min{ ggg^5g , fe^}. Thus (£,£) T (£) 
is integrable at -co. 

Step 2. We show that the function iwa(C)) =: J*ioo(^( s )> v(s))ds is 

integrable for £ close to -co. By step O (wi, w 2 )(Q is well defined and 
convergent. 

Let 5 = max _ {u + k\ — a\v, r(a 2 u + k 2 — v)} + 1 and (3 = 5/c. 

0<u<aik2,0<v<k2 

Firstly we show the function 

u 



, , (0 



is monotonically decreasing. 
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Since 



*~ K ( ~ ) (0 = e-« ( - - £ ) (0 



?! 



t> 



1 -p£ ( —J*u + u + (—ki—u){aiv — u) — 5u 



c I — J * v + v + r(^2 — v){— a 2 u + *y) — Sv ' ' s ' 



< 0, 



it follows that the function e ^ f ~ ^ (£) is decreasing. 



■it 



v 



It is easy to see that J * I ~ (s)ds is also integrable on (— oo,£) for £ 



sufficiently negatively large. Integrating (|3.7[) from — oo to ^ we get 



k\ — e —a\k\ — £ \ I ( u 
-a 2 rk 2 - e rk 2 - e J J ^ \ v 



(s)ds 



Integrating this inequality again from y to £ we get 



k\ — e —a\k\ - e \ n ( w\ 
-a 2 rk 2 - e rk 2 - e J J y \ w 2 



(s) 



^jf(l)^ + jf( J -(S)-(Sll* 



Since J has compact support, we can again use Fubini's Theorem and 
Lebesgue's Theorem to obtain 



as y — t- — oo. Thus ( ^ J ( s ) * s integrable on (— oo,£). 

Step 3. The function (u,v) decays at least exponentially at — oo. 
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Similar to step 2, we have 



k\ — e —a\k\ — e \ f ( w\ 



-a 2 rk 2 - e rk 2 - e J J_ OD \ w 2 



< c 



Wl 

w 2 

-),o + /; J(W( )((-) K ) + f' (»)«*)* 

< (c+ fj\J(t)dt) ( ^ )(«) + £( I) K) fj\J(t)(e-"-\)dL 

Here we have used the facts that ( ~ ^ (£) > and e _/3 ^ ^ ~ ^ (£) is 

decreasing. 
Furthermore, 







/3 V v J J-oc V -«2rA;2 - e rk 2 - e J J_ x ^ w 2 7 

< ( c + £ mm) ( £ ) (0 + 1( f ) «> /I 

Hence there exists a positive matrix M (see step [3]) with all positive entries 
such that 

/1((S) W+ (") W )*^<(S) K)+ (") K) ) 

Noting that ( ~ ) (s) and ( Wl ) (s) are increasing on (— oo,£), we have 



v I \ w 2 



Thus for large r, there is some < k < 1 such that 



?/>! 
W 2 



(f-r) + (»)tt-r) S *((S)(0 + («) (e )) 
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Let H(() = e-* ( ( ^ ) (f > + ( ~ ) («)) . where 1 = 7 '» i 

HK - r) = e -7tt-)((»i) (? _ r)+ (5) (f _ r) ) 

< ^((^ )«)+(?)( 0)-*K) 

The fact that < u < aiA;2 and < u < k 2 implies 



W 2 



and thus lim = 0, and H(£) is bounded. Consequently we have 

Zl )(£)+( |)(0<O(e^)as^-oc 
and therefore 

' ~ ) (0 < 0{e«) 

for £ — >• — oo. □ 

The inequality (u,v)(£) < 0(e yx ) says that the two components of the 
traveling front may approach the equilibrium with different rates, however, 
we will show this will not happen and give the exact decay rates of the front 
solutions. Recalling (see the proof of the previous Lemma) that at — oo 
the limit matrix of the Jocobian of the vector functions ((— k\ — u)(a,\v — 
u),r(k 2 -v)(-a 2 u + v)) is 

(k\ —a±ki 
-a 2 rk 2 rk 2 
and it has two positive eigenvalues 



ki + rk 2 ± yj (k\ + rk 2 ) 2 - <irkik 2 (l - a\a 2 ) 
Mi,2 = 2 > U - 

For any c > each of the following equations 

M n = / J ( s ) e Xs ds+c\- l-^i = 0, M 22 = / J(s)e Xs ds+c\-l-fi 2 = 
Jr Jr 

has exactly one positive root. For c > c* (c = c*) let /2(c) (/2(c*)) be the 

smaller positive root of the above equations. We also denote the smaller 

positive root of the equation 

M_(A)= f J{£,)e^dt-c\ + r{l-a 2 )-l 
Jr 
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by Ai(c) for c > c* and its double positive root as A* for c = c*. 

Theorem 12. For every c > c* , the traveling wave solution as derived in 
Theorem^ has the following asymptotic properties: 
(i). Corresponding to the wave speed c > c* , 



as £ — > — oo; and 



(3-9) = r - T e-^+o(e 



12 



as £ — >• +oo, where A\, A2, Ai, A2 are positive constants, 
(ii). Corresponding to the wave speed c = c* , 

(3.10) f:«n=f^t-M 



v(0 J \ A21 + A 22 i 
as £ — > —00; and 

(3 ' u> ("IMt 2 )- (I: 

as £ — > +00, where A12, A 2 2 < 0, An, A21 G M and An, A21 > 0. 

Proof. For every c > c*, the traveling front solution (u(£,),v(£,)) T to the 
equation (jl.5p satisfies 

where (u,v) and (u, are lower and upper solutions of (|1.5p . 

Lemma implies that the upper- and lower-solutions as derived in Lemma 
[3]and Lemma H] have the same asymptotic rate at £ — > —00. Then (|3.8[) and 

(pnop follow. 

We next prove (|3.9p and (|3.1ip . By Lemma [T0| the traveling front ap- 
proaches the equilibrium {a\k2, fo) at least exponentially, now we determine 
the exponential rate. It is convenient to start with the derivative of the 
traveling wave (wi,W2)(0 = ( u \ v> )(0 which satisfies the system 

{J * v\ — v\ — cv[ + (2u — a\v — \)v\ = a\(u — l)t>2, 
J * V2 — v 2 — cv' 2 + [r(l — 02) + ra2U — 2rv]v2 = —a2rvv\ 

It is easy to see that v\{£) can not decay slower than t>2(£)- For, if 
we divide both sides of the second equation of (13.12p . and take limit for 
£ — > —00, we would have a finite number in the left hand side and infinity 
at the right hand side. A contradiction. Similarly we can show that V2 can 
not decay slower than v\ by the first equation of (|3.12p . This implies that 
a\k,2 — u = 0(/c2 — v ) as £ — > +00. 
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We can change the matrix M into diagonal by Jordan canonical transform, 
in fact let 



P 

then 



— a\k\ —a\k\ 
Ml ~ k i 1*2- k i 



i / M2 - ki ki - //I \ / [ii 

P- 1 = - so that P- 1 MP = 

aihfa-K) \ aikl _ aih ) \ o 

Now writing system (|3.4|) in vector form with U = (u, v) and F = ((— k\ — 
u){a\v — u),r(k2 — v)(—a2U + v)): 

J*U -U + cU' - F{U) = 0. 
Introduce the transformation 







f wi \ 


[ V I 










\W 2 



(3-13) I j(f) = pl j(^), 

then W = {w\,W2) satisfies the following equation 

J *W -W + cW' - P- X F(PW) = 

or equivalently 
(3.14) 

/ m o \ ( mi o \ 

J *W -W + cW - \W = P- X F{PW) - I W. 

V fl 2 J \ fx 2 J 

It is easy to verify that as £ — > — oo, 

/Mi 0\ 

(3.15) P _1 F(PW0- \W = 0{w\ + wl) 

\ M2 / 

Since w\ and ^2 are linear combinations of u and we have (wi,W2) < 
0(e 7 ^) as £ — > — 00. It then follows that for < ReX < 7 the two side 
Laplace transforms 

Wx(X) = [ e- x twi(t)(% and W 2 (X) = [ e^ii^H 
are well defined. Note that 
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J*W(£)d£ = [ J(s)e Xs f W(C + s)e' x(ii+s) dCds 
Jr Jr 

Wi 



w 2 



(A) / J(s)e As ds 



we then have 



/ W 1 (X) \ . ( Ml . 

(3.16) M(c,A) = / e-^(P^F(Plf)- W)d£ 

\ W 2 {X) I Jr \ M2 / 



M(c, A) 



where the matrix M is given by 

J K J(s)e Xs ds + c\-l- m 

J R J(s)e Xs ds + cA - 1 - /i 2 

and from f)3. 15|) we see the right hand side of (|3.16p is well defined for 
< ReX < 27. According to [p 58 Theorem 5a of D. V. Widder], there is 
a positive A such that Wi(A) and ^(A) are analytic for < ReX < X and 
one of W\ and W 2 has singularity at A = A. This shows that for c > c*, W\ 
and W 2 are defined for ReA < p, = min{/2i , fi 2 } where fii , i = 1,2 are the 
positive root of the two equation: 

(3.17) Mn = / J(s)e Xs ds + cA - 1 - m = 0, 

it 

and 

(3.18) M 22 (A) = / J(s)e Xs ds + cA - 1 - /i 2 = 

Jr 

respectively. 

Now Let P- 1 F(PW) = (F^w^wz), F 2 (w-l,w 2 )) and write 
(3.19) 

wi(-0e ? df = — . / Wi(Qe ^d£, 

Mn(c, A) Jo 

and 
(3.20) 

[ + °° t c\ Km ke- x ^[F 2 ( Wl ,w 2 ) - ^ 2 w 2 ]dj f°° _ A . 

/ ^2(-0 e = — : rr / ™ 2 (£)e 5 d6 

Jo M 22 (c, A) Jo 
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Observing that the second terms wi^(O e X ^d£, on the right hand side 
of (|3.19p and (|3,20p are analytic for ReX > 0, and there is no positive root 
of (ETT71) ( (pTTHjl ) with ReX = fix (ReX = /2 2 ) other than A = Mi (A = fi 2 , 
respectively). To see this ([H [H]), we let A = Ml + ifi, /3 €E R such that 



= Mnifix+iP) 

= I R e^+Wyj(y)dy - c(/2i + ip) - 1 - Mi 

= /r eAiy J (y) cos - c/Zi - 1 - Ml + (/ K e^ lS/ J 0) sin /3ydy - c/3)i 

= / R e^y.J(y){-2 sin 2 f + l)dy - cfix - 1 - Ml + (/ R e^J(y) sin/?ydy 

= -2 J R e*»J(y) sin 2 f + (/ K e^J(y) sin /3ydy - c/3)i 



therefore f3 = 0. We can also show that M22(A) has a similar property. 

Hence we can apply the modified Ikehara's Tauberian Theorem (see [I]) 
to the function wiand W2 separately to get 

Wi = A ie M + o(e^), i = 1,2 

for Ai>0,i = 1,2. 

By (|3.13p , u and v are linear combinations of w\ and W2 ■ We see that they 
both decay exponentially at a rate of p, = min{//i,M2} at — oo. Switching 
back to u, v from u and v and reversing back the sign of £, we have the 
conclusion of the theorem. □ 

We introduce a comparison principle for non-local system that will be 
used in showing the uniqueness of the traveling front solution of (II. 5h . 

Lemma 13. Let the C 1 vector function U(£) = (tti(£),«2(£), ' ' ' >^n(£)) T 
and U_(^) = («i(£)>M2(£)' " " " jl*n(0) T ^ e monotonically increasing in M anc? 
satisfy the following inequalities 



which implies that 




D(Z7) - cU' + F(C7) < D(U) - cU' + F(U) for£ £ [-N, N] 



and 



U(-N)<U(0, U(0<U(N) forte(-N,N] 



U(0 < U(0 for £ 6 (-oo, -N] U [N, +oo), 
where D(U) is a diagonal matrix with entries Di(Ji*Ui—Ui), Di > and J{ is 
a positive, even integration kernel with unite mass fori = 1, 2, • • • , n,F(U) = 
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(Fi(U), ■ ■ ■ ,F n (U)) T is C 1 with respect to its components and |^ > for 
i + 3\ h 3 = 1,2,- •• ,n, then 

LL(0<U(0 fort€[-N,N]. 
Proof. See [16\ Lemma 12]. □ 

Theorem 14. Assume that < a±,a2 < 1 and (|2.7p . for every c > c* , 
system (jl.5p /ias corresponding a unique (up to a translation of the origin) 
traveling wave solution. 

Proof. We only prove the conclusion for traveling wave solutions with as- 
ymptotic rates (|3.8|) and (j3.9|) . since the other case can be dealt similarly. 
Let Ui(£) = M£),vi(£)) T and U 2 (£) = (MOMOf be two traveling 
front solutions of system (|1.5p with the same speed c > c*. Then there exist 
positive constants Ai, Bi,i = 1, 2, 3, 4 and a large number N > such that 
for $ < —N, 

(3-21) ^^(te^)^^ 

(3-22) C/2 ^=(teS) + ° (eAl?); 
and 

(3-23) ^K)=(t-"^« ) +0(eA * £) 



(3.24) um = ( °*_- B ^f ) + 

The traveling front solutions of system (jl.5p are translation invariant, thus 
for any 9 > 0, f/f (£) := + 0) is also a traveling wave solution of (jl.5p . 
By (|3.2ip and (|3.23p . the solution f7i(£ + 0) has the asymptotic rates 

a ( Aie Xld e Xll i \ i * 

(3-25) ^)=(t^V 1 + ° ( } 

for £ < —N; 

for | > AT. 

Choosing > large enough such that 
(3.27) ^ie Aie > A 3 , 



(3.28) 



THE COMPETITION INCLUSION 



21 



(3.29) B ie X:ie < B 3 , 

(3.30) B 2 e X3d < B 4 . 
then one has 

(3-31) C/f(0 > U 2 (0 

for £ G (-oo, -iV] U [N, +oo). 

While on the interval [— JV, N] we can apply Lemma[T3lto have Iff (£) > ^(C) 
for £ G [-iV, iV]. This shows that 

C/f(C)>C/ 2 (£), (GK. 
Now, decrease # until one of the following situations happens. 

(i) There exists a 9 > such that C/f (£) = ^(O- in this case we have 
finished the proof. 

(ii) There exists a 6 > and £i G M such that one of the components of 
U e and £/ 2 are equal there; and for all £ G M, we have Uf (£) > t^2(£)- One 
applying the Maximum Principle on R for that component, we find Iff (£) 
and U2{C) must be identical on that component. Without loss of generality, 
we suppose that the component is the first component. Then Uf (£) — ^(C) 
satisfies ()3. 1|) and (j3.3|) . Plugging w\ = into (|3.ip again we find that 
there is at least one £^ such that ^2(67) = 0. Then by applying Maximum 
Principle to (|3.2[) . we have tU2(£) = for ( £ 1. We have then return to 
case (i). 

Consequently, in either situation, there exists a 9 > such that 

for all £ G M. □ 

The next theorem shows that lower bound c* for the wave speed c is 
optimal, hence it is the critical minimal wave speed. 

Theorem 15. There is no monotone traveling wave solution of (|1.5|) for 

any < c < c* . 

Proof. The proof is basically the same as that of 0] so we skip it. □ 
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